Spin-Hall interface resistance in terms of Landauer type spin dipoles 
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We considered the nonequlibrium spin dipoles induced around spin independent elastic scatterers 
by the intrinsic spin-Hall effect associated with the Rashba spin-orbit coupling. The normal to 
2DEG spin polarization has been calculated in the diffusion range around the scatterer. We found 
that although around each impurity this polarization is finite, the corresponding macroscopic spin 
density, obtained via averaging of individual spin dipole distributions over impurity positions is 
zero in the bulk. At the same time, the spin density is finite near the boundary of 2DEG, except 
for a special case of a hard wall boundary, when it turns to 0. The boundary value of the spin 
polarization can be associated with the interface spin-Hall resistance determining the additional 
energy dissipation due to spin accumulation. 
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I. INTRODUCTION 

Most of the theoretical studies on the spin-Hall effect 
(SHE) has been devoted to calculation of the spin cur- 
rent (for a review see Ref. [l|). Such a current is a linear 
response to the external electric field E which induces 
a spin flux of electrons or holes flowing in the direction 
perpendicular to E. This spin flux can be due either to 
the intrinsic spin-orbit interaction (SOI) inherent to a 
crystalline solid 2 , or to spin dependent scattering from 
impurities^. The spin-Hall current, as a response to the 
electric field, is characterized by the spin-Hall conduc- 
tivity. On the other hand, similar to the conventional 
Hall effect, one can introduce the spin-Hall resistivity 
from calculating the local chemical potential difference 
jj, s = jj, t — A* I in a response to the DC electric cur- 
rent. For 2DEG in a local equilibrium this potential 
difference can be related to the z-component (perpen- 
dicular to 2DEG) of the spin polarization, according to 
S z = Npfi s , where Np is the density of states near the 
Fermi level. Therefore, the spin-Hall resistivity is closely 
associated with spin accumulation near interfaces. It 
should be noted that measuring spin polarization is thus 
far the only realistic way to detect SHE^i 5 -. For inter- 
faces of various nature such an accumulation has been 
calculated in a number of works^ 7 ^ 9 ' 10 ? 11 ! 12 . A typical 
example to study spin accumulation is an infinite along 
x-direction 2D strip with a width w along y-direction. In 
this geometry the DC current flows in x-direction, while 
the spin-Hall current flows in y-direction with the spin 
density accumulating near boundaries. An analog of the 
Hall voltage could be a difference of fi s on both sides of 
the strip. There is, however, a fundamental distinction 
from the charge Hall effect. In the latter case, due to the 
long-range nature of the electric potential created by con- 
serving electric charges, the Hall voltage is proportional 
to the width of the strip. In contrast, the spin-Hall elec- 
trochemical potential at the interface does not depend 
on w as w — ► oo because spin relaxation essentially sup- 
presses the long range contribution to spin polarization 
buildup near interfaces. Hence, it is sensible to introduce 



an interface spin-Hall resistance, which is the proportion- 
ality coefficient between the interface value of /i s and the 
electric current density. 

Below, we will consider the spin-Hall resistance from 
the microscopic point of view. This approach is based 
on the Landauer'si^ idea that at a given electric current 
each impurity is surrounded by a nonequlibrium charge 
cloud forming a dipole. Combined together these dipoles 
create a voltage drop across the sample. Therefore, each 
impurity plays a role of an elementary resistor. In a sim- 
ilar way, nonequlibrium spin dipoles could be induced 
subsequent to the spin-Hall current. One may expect 
that the spin cloud will appear around a spin-orbit scat- 
terer in case of extrinsic SHE, as well as around a spin- 
independent scatterer, in case of the intrinsic effect. The 
latter possibility for a 2D electron gas with Rashba inter- 
action has been considered in Ref. 1 14 . The perpendicular 
to 2DEG polarization was calculated within the ballistic 
range around a scatterer. On the other hand, in order to 
study spin accumulation and the spin-Hall resistance on 
a macroscopic scale, one needs to calculate the spin den- 
sity distribution at distances much larger than the mean 
free path / of electrons. Below, we will extend the Green 
function method of Ref. to the diffusive range. In 
Section II the spin density distribution around an indi- 
vidual target impurity will be calculated. In Section III 
we will consider the interface spin accumulation created 
by spin dipoles randomly but homogeneously distributed 
in space. A relation between spin-Hall resistance and en- 
ergy dissipation will be discussed in Sec. IV. A summary 
and discussion of results will be presented in Section V. 



II. SPIN CLOUD INDUCED BY A SINGLE 
IMPURITY 

As known, the electric field applied to a homogeneous 
2DEG with Rashba SOI induces a parallel to 2DEG com- 
ponent of the nonequlibrium spin polarization 15 . The 
spin-Hall effect produces, however, a zero spin polariza- 
tion in its z-component. This understanding about such 
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a homogeneous gas has implied an averaging over impu- 
rity positions. An impure system, on the other hand, 
can not be uniform on a microscopic scale. The effect of 
each impurity on the spin polarization could be singled 
out by considering an impurity (a target impurity) at a 
fixed position while taking at the same time the aver- 
age over positions of other impurities. In such a way the 
Landauer electric dipole has been calculated 16 ' 17 . The 
electron density around a target impurity represented 
by the elastic scatterer was found from the asymptotic 
expansion of the scattered wave functions of electrons. 
At the same time, wavevectors of incident particles were 
weighted with the noncqulibrium part of the Boltzmann 
distribution function. We will employ another method 
based on the Green function formalism 14 . Within this 
method the spin density response to the electric field E 
is given by the standard Kubo formula with the scatter- 
ing potential of the target impurity incorporated into the 
retarded and advanced Green functions G r / a (r, r', w) de- 
noted by the superscripts r and a, respectively. As such, 
the n-component of the stationary spin polarization is 
given by 



dV 



du> drip (to) 
2n dw 



Tr[a n G r {v, r', w)(vE)G«(r', r, w)] , (1) 

where the overline denotes averaging over impurity po- 
sitions, the trace runs through the spin variables and 



np{uj) is the Fermi distribution function. To avoid fur- 
ther confusion we note that the angular moment is ob- 
tained by multiplying S n {r) by h/2 and e is the particle 
charge, which is negative for electrons. At low temper- 
atures only uj in close vicinity around Ep contributes 
to the integral in (fTJ). Therefore, below we set u = Ep 
and omit the frequency argument in the Green functions. 
Further, v is the particle velocity operator containing a 
spin dependent part associated with SOI. Writing SOI in 
the form 



H so — hk • cr 



one obtains the velocity operator 



(2) 



m 



9hk • (r 



(3) 



where <r=(cr x , a y ,a z ) is the Pauli matrix vector. In case 
of Rashba interaction the spin-orbit field hk is given by 



h x = aky , h v = -ak x . 



(4) 



We assume that the target impurity, located at r^, is rep- 
resented by a scattering potential U(r — rj). The Green 
functions in ([T]) have to be expanded in terms of this 
potential. Up to the second order in U one obtains 



G r/a (r,r') = G r/a(0) (r,r') + / ds 2 G r/a(0) (r, s)U(s - n)G r ^(s, r') + 

J ds 2 ds' 2 G r/a{0) {Y,s)U{s - r i )G r/a(0) (s,s')E/(s / - r % )G r/a{Q) (s' ,r') . (5) 

I 



The unperturbed functions G r / a ^ depend, nevertheless, 
on scattering from background random impurities. The 
latter create the random potential V sc (r) which is as- 
sumed to be delta correlated, so that the pair correlator 
{V sc (r)V sc {r')) = T5(r - r')/nN F , where T = l/2r is ex- 
pressed via the mean elastic scattering time r. The delta 
correlation means that the corresponding impurity po- 
tential is the short range one. In fact, the potential of the 
target impurity could be different from that of the ran- 
dom impurities. It might be a special sort of impurities 
added to the system. On the other hand, the target and 
the random impurities would be identical if one would 
try to employ the spin dipoles for the interpretation of 
spin accumulation near interfaces. 

After substitution of ([5]) into (fTJ) one must calculate 
background impurity configurational averages contain- 
ing products of several Green functions G^ ' . Assuming 
that the semiclassical approximation Epr 3> 1 is valid, 



the standard perturbation theor y 18 ' 19 can be employed 
whose building blocks are the so called ladder perturba- 
tion series expressed in terms of the unperturbed average 
Green functions 

Gl/ a = J d 2 (r - r')e ik(r - r,) G r /°( )(r, r') (6) 
given by the 2x2 matrix 

G r k /a = (E F -E k -h k -<r± iry 1 , (7) 

where E k =k 2 / (2m*) . When averaging the Green func- 
tion products, within the ladder approximation only 
pairs of retarded and advanced functions carrying close 
enough momenta should be chosen to become elements 
of the ladder series. After decoupling the mean products 
of Green functions into the ladder series, the Fourier ex- 
pansion of (TT]) can be represented by diagrams shown 
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FIG. 1: Examples of diagrams for the spin density. Scattering 
of electrons by a target impurity is shown by the solid circles. 
Dashed lines denote the ladder series of particle scattering by 
the random potential, p, k, k' are electron momenta. 

at Fig.l. In these diagrams the diffusion ladder renor- 
malizes both the lefthand and righthand vertices. The 
renormalized lefthand vertex E z (q) is associated with the 
q-th Fourier component of the induced spin density, and 
the corresponding diffusion propagator enters with the 
wavevector q. In its turn, the righthand vertex T(p) 
related to the homogeneous electric field, is represented 
by the ladder at the zeroth wavector. The correspond- 
ing physical process is the D'yakonov-Perel spin relax- 
ation of a uniform spin distribution. This lefthand ver- 
tex alone contributed to the ballistic case result, while 
S 2 (q) has been taken unrenormalized due to large values 
of q ^ 1/(uft)1 in the ballistic regime. Figs. 1(e), 1(f) 
represent some diagrams where the diffusion process sep- 
arates two scattering events. As it will be shown below, 
such diagrams give rise to small corrections to the spin 
density, and can be neglected. Hence, the main contribu- 
tion comes from diagrams similar to those in Figs. 1(a)- 
1(d). The corresponding spin polarization has the form 

^(<l) = ^ E 7MG£ k £ 2 (q)G k+q , p T(p)] . (8) 
P ,k 

The functions G k ( k are formally represented by the 
Fourier expansion of Eq. ([5]) with respect to r and r', 
providing that the respective average values G(°)(r, r') 
are substituted instead of G^ (r, r'). Evaluating the pair 
products of such functions in (jSJ , one should take into ac- 
count only terms up to the second order with respect to 
the scattering potential U. 

The vertices S z (q) and T(p) can be easily calculated. 
As it was discussed in Ref. LU], due to considerable can- 
cellation of diagrams which is known from literature on 
the spin-Hall effect, T(p) acquires a quite simple form in 
a special case of Rashba SOI. Namely, 

T(p) = -^p-E. (9) 



In its turn, £ z (q) is expressed in terms of the diffusion 
propagator. Indeed, let us represent this vertex using 
a basis of four 2x2 matrices r° = 1 and t % — a 1 with 
i = x, y, z. Then, £ z (q) can be written as 

£ z (q) = ^ zb (q)T b , b = 0,X,y,Z (10) 
b 

where D zb (q) are the matrix elements of the diffusion 
propagator satisfying the spin diffusion equation, as it 
was described in Ref. and references therein. The non- 
diagonal element D z0 (q) appears due to the spin-charge 
mixing and it is zero for SOI of quite general form, in- 
cluding the Rashba interactio n 20 ! 22 ^.. Finally, from (jHJl, 
using ([9]) and (fTO]) , we express S z (q) in the form 

S z (q)= J2 D* n (q)I n (q), (11) 

n—x,y,z 

where / n (q) 

/n(q) = 2^ 5> ' E )^[^ k ^G k+q ,p] ■ (12) 
P .k 

The function /™(q) has a simple physical meaning. For 
n = x,y,z it represents a source of spin polarized par- 
ticles emitting from the target impurity. Their further 
diffusion and spin relaxation result in the observable po- 
larization. This source term feature is conceptually simi- 
lar, though different in its context, to the original charge 
cloud consideration when SOI is not present and Boltz- 
mann equation is used to describe the subsequent back- 
ground scatterin g 13 ! 16 . For q <C l^ 1 <C &f the source can 
be expanded in powers of q. Therefore, the wavevector 
independent terms represent the delta source located at 
Ti, while the linear in q terms are associated with the 
gradient of the delta-function. Below we will keep only 
the constant and linear terms for each n-th component 
I n (q) and assume, for simplicity, the short range scatter- 
ing potential f/(r), so that its k-th Fourier transform is 
simply J7exp(— ik • i - ;), where U is a constant. Further, 
I n (q) can be written as 

/"(qHW+W, (13) 

where I\ and I2 are of the first and the second order 
with respect to the scattering potential U, respectively. 
Accordingly, I x and I2 are represented by Figs. 1(a), 1(b) 
and Figs. 1(c), 1(d), respectively. Using (O to express 
Green functions G k ( k in (JT2J) we obtain 

/l(q) = ^"'"E^x 
p 

Tr[G r p G a p (* n G' p+cl + Gl_^ n )] , (14) 

and 

pk 

(G£a"G k+q - 7 ^ +q + 7 G»_ q a")] , (15) 
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where 



7 = ilm[) G£] = iirN F 



k 



(16) 



In our following consideration we let the x-axis to be par- 
allel with the electric field, and the z-axis to be perpendic- 
ular to the 2DEG. The system Hamiltonian is symmetric 
under a symmetry operation combining a reflection from 
the plane perpendicular to the y-axis, that means p y — > 
— p y , and a unitary transformation a' 1 — > a y a l a y . Ap- 
plying this transformation to (fT2|) one can easily see that 

r x (qx,q v ) = ~P(q x ,~q v ), i z (q x ,q v ) = -i z (qx,-q y ) and 

I y (q x ,qy) — I y (qx,—qy)- Making use of another symme- 
try operation p x — > p y — > -p y and er 2 — > cr z crV z , 
we obtain I x (q x ,q y ) = P(-q x ,-q y ), P(q x ,q y ) = 
-P(-q x ,—q y ) and P(q x ,q y ) = P(—q x ,—q y ). From 
these relations it is easy to see that expansion of P into 
power series starts from linear in q terms, while the lead- 
ing term in I y is const and the next one is quadratic in 
q. On this reason only const will be taken into account 
in I v . The expansion of P starts from q x q y , and this 
source component will be neglected. 

Calculation of Ii and I 2 given by Eqs. (fl~4")) , (fT5j) is 
based on the standard linearization near the Fermi level, 
thus ignoring band effects giving rise to small correc- 
tions ~ hk F I Ep,T/ Ep. Further, the diffusion approx- 
imation is valid at q -C l/l. At the same time, the 
characteristic lengthscale is determined by the spin re- 
laxation length l so , which is the distance a particle dif- 
fuses during the D'yakonov-Perel' spin relaxation time 
t so = 4(/i 2 . f t) _1 . The corresponding diffusion length 

lso = V Dt so , where D = v f t/2 is the diffusion constant. 
Hence, l so = vf /hk F - Taking q ~ l/l so one finds that the 
diffusion approximation is valid if h^ F /T <C 1. Therefore, 
within this approximation we will retain only the leading 
powers of hk F /T < 1. In such a way, direct calculation 
of 7f with the Green functions and SOI given by Eqs. 
([7|) and (j4|), respectively, shows that both 1\ and If are 
small by a factor T/Ep. For example, using the relation 



9 n r/a 

dEp Uk 



(17) 



which follows from ([7]), evaluating If at q — 0, one can 
represent the corresponding sum in (|14[) as 



i) 



8E f 



Y,PxTr[G r p G a p a y ] = 



-dE-p {- NFm dp-) 



(18) 



linear in U contribution to the spin density was associ- 
ated with fast Friedel oscillations. It is clear that their 
Fourier transform will be small in the range of q <C kp 1 . 



At the same time I y and 7| 



given by 



P = 

P = 
P = 



v d Npm*ah 2 kF ^ 

V 

-iq yVd N F hl F — 
0, 



are not zero. They are 



r 



(19) 



where V = irNpU 2 and Vd — eEr/m* is the electron 
drift velocity. If the target impurity is represented by 
one of the random scatterers, we get T' = T/ni, where 
Ui is the density of impurities. 

In the above calculation we did not take into account 
the diagrams shown in Figs. l(e)-l(f) and those similar 
to them. It can be easily seen that such diagrams contain 
7™ as a factor. For example, the sum of diagrams at Fig. 
l(e)-l(f) contains as a multiplier the sum of diagrams 
shown at Figs. l(a)-l(b). Therefore, such diagrams are 
small by the same reason as J™ are, at least, in the most 
important range of / <C P 1 ■ Particularly in this range 
of small / the diffusion propagator between the two scat- 
tering events in Figs. l(e)-l(f) becomes large. 

Now one can combine the source P with the diffusion 
propagator to find from Eq. pip the shape of the spin 
cloud around a single scatterer. Taking into account (|19p . 
Eq. (fTTj) is transformed into 

S z (q) - -v d N F hl F ^ HjD^iq) - 2mW(q)] . 

(20) 

The matrix elements D y (q) satisfy the spin diffusion 
equatio™^! 



53 [ -5 a Dq 2 - Y a + i R am qm J P> 13 (q) - -2F^ , 

l \ rn J 

where the matrix P l determining the D'yakonov-Perel' 
spin relaxation rates is given by 



^(5 u hl F ~P kp h l k 



(22) 



with angular brackets denoting averaging over the Fermi 
surface. In the case of Rashba SOI, Eq. ((4]), one gets 
T zz = 4rh 2 k and T xx = T vy = 2rh\ . The last term in 



lhs of (|2Tj) is associated with spin precession in SOI field. 
It has the form 



R ilm =4r ^ £ «P ( h P t 



(23) 



In case of Rashba SOI with the constant coupling 
strength a and energy independent parameters T, m* and 
Np, the sum (fT5]) is equal to 0. Otherwise, it is finite, 
but small due to the smooth energy dependence of these 
parameters. Similar analysis, although not so straight- 
forward, can be applied to If , which is linear in q. The 
smallness of P can be also seen from Ref. fLJ where the 



For the Rashba SOI the nonzero components are 



AiDm*a qi . (24) 



We ignored in (|2"Tj) a small term which gives rise to the 
spin-charge mixing£ i 22 i 23 . This mixing is already taken 
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III. SPIN ACCUMULATION IN A 
SEMIINFINITE SYSTEM 




FIG. 2: Spatial distribution of S z component of the spin den- 
sity around a single scatterer. The unit of length = l so 



into account in the source term because I n for n — x,y, z 
describes the source of the spin polarization in response 
to the electric field. From Eqs. (|2T]) -(j24 |) one finds 



D zz 
-D zy 
jyvv 



<f + l 



2/4 f t 2 {q 2 + 2){q 2 + 1) - Aq 2 
1 2iq y 



D v 



1 



2h 2 p T 2 (q 2 + 2){q 2 + l)-Aq 2 
q 2 + 2 



2h 2 T 2 (q 2 + 2)(q 2 + l)-Aq 



in? 



(25) 



In this section we will consider a semiinfmite electron 
gas y > bounded at y = by a boundary parallel to the 
electric field. Our goal is to calculate a combined effect 
of spin clouds from random impurities. It is important to 
note that summation of spin dipoles from many scatter- 
ers does not result in a magnetic potential gradient in the 
bulk of the sample. This is principally different from the 
Landauer charge dipoles which are associated with the 
macroscopic electric field. The origin of such a distinc- 
tion can be immediately seen from (|26[) . The magnetic 
potential, as it was defined in Sec. I, is proportional to S z . 
Taking its gradient one gets q y S z . After averaging over 
impurity positions q — > 0, so that q y S z — ► 0. It happens 
due to spin relaxation, which provides at q = a finite 
value of the denominator in (|26p . At the same time, in 
case of the charge cloud, the denominator of the particle 
diffusion propagator is proportional to q 2 . Hence, the 
corresponding gradient of the electrochemical potential 
(electric field) is finite at q = 0. Although the bulk mag- 
netic potential is zero, one can not expect that it will also 
be zero near an interface. In order to calculate the spin 
polarization near the boundary, Eq. (|2ip with q = — iV 
and 2TS(r)5ij in the r.h.s. has to be solved using ap- 
propriate boundary conditions. With the so obtained 
D^(r), the resultant spin density induced by impurities 
placed at points is given by Eq. lTTj) 



c')I? ot (r'), (28) 



where 2q = l so q denotes the dimensionless wavector. 
Substituting (55]) into (201) we finally find 



S z = 



m*a r' 
-2iv d ——N F - 



<7 y (<? 2 +3) 



r (q 2 + 2)(q 2 + l)-4q 



z2 



(26) 



and 



m*a V 
S y = 2v,i^—Nf 



(iq 2 + 2) 



(27) 



r (q 2 + 2){q 2 + 1) - 4q 2 ' 

To restore the conventional units we added h into 
(f27|) . The z component of the spin density in real space 
is shown in Fig. 2. According to expectations it has the 
shape of a dipole oriented perpendicular to the electric 
field. Its spatial behavior is determined by the single 
parameter l so , which gives the range of exponential de- 
cay of the spin polarization with increasing distance from 
an impurity. The S y component averaged over impu- 
rity positions gives the uniform bulk polarization. It is 
interesting to note that when the target impurities are 
identical to background ones (r' = T), the so obtained 
uniform polarization S y \ q ^o coincides with the electric 
spin orientation^ 5 S y = 2vdm*aNp /h. 



where the source term is obtained by the inverse Fourier 
transform of (fT9l) : 



^L( r ) = v dN F m*ah 2 kF -^—^2S( 

1 Tli 



(29) 



where the relation r" = T/m is used because we assumed 
that the target impurities are identical to the random 
ones. The macroscopic polarization is obtained by av- 
eraging of (j28f and ([29]) over impurity positions. After 
averaging over Xj, and the semiinfmite region y^ > the 
spin polarization source (|29|) transforms to I£ v (y): 



!L(y) = v d N F m*ahl 



1 



-v d N F h 2 k J(y-0+) 



2T 2 



(30) 



It follows from (f28|) that the corresponding mean value 
of the spin polarization S av (y) satisfies the diffusion equa- 
tion (f2"Tj) with the source 2r/™ v (y) in its rhs. The so ob- 
tained diffusion equation, however, is not complete. One 
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should take into account that the boundary itself can 
create the interface spin polarization. Most easily it can 
be done in the framework of the Boltzmann approach. 
In terms of the Boltzmann function the spin density is 
defined as S av (y) = X)k§k, an d the charge density as 
Skffk- The equation for the Boltzmann function can be 
written in the form (see e.g. Ref. [25h 



diffusion equation which acquires the form 



VyVygk + 2(g k x h k ) + eE- 
= - (S E (y) - gk) , 

T 



(31) 



(0) 



where S E (y) = 5{E - E F )S av (y)/N F and g^ 
— hkS(E — Ef) is the equilibrium Boltzmann function. 
The terms proportional to the charge component of the 
Boltzmann function have been omitted in pip due to the 
system local electroneutrality, at least in the scale of the 
mean free path, which is the smallest characteristic scale 
of gk spatial variations. The spin polarization source as- 
sociated with the boundary is given by a direct action of 
the electric field, without taking into account secondary 
scattering from impurities. Hence, the term with S av (y) 
in r.h.s. of Eq. (|3"Tj) can be ignored. Also, the boundary 
independent bulk part of gk has to be subtracted from 
the general solution of pip . The so obtained interface 
Boltzmann function will be denoted as gkif- The corre- 
sponding spin density is Sif(y) = X)k§kif- I n order to 
calculate gkif, Eq. (|3Tj) has to be supplemented with the 
boundary condition. For a hard wall specularly reflecting 
boundary the condition is simply 



=0 — gkx,-k v \z=0 ■ 



(32) 



This condition means that the spin orientation does not 
change after specular reflection from the interface. The 
solution of Eq. (|3"Tj) satisfying (f3"2")) can be easily found. 
As a result, up to o(a 2 ) we obtain 

SfM = S§(y) = 

Sg(y) = 8v d a 2 rm* M(^k - E F )e~ v & . (33) 

ky>0 

Within the diffusion approximation the second of these 
equations represents a delta source of the spin polariza- 
tion with intensity 



1 



1 

dyS? i (y)=v d N F h 2 kF -. (34) 



This source is exactly of the same magnitude, but oppo- 
site in sign to the spin polarization emerging from im- 
purities, that is represented by the integral of 2TI* v (y), 
with I^viy) given by Eq. ([30)1 . Taking into account that 
both sources are located at the interface, so that they 
cancel each other, one sees that only y-component of the 
source originating from impurity scattering retains in the 



dy 2 
d 2 S y 
dy 2 



Am* a 



+ Am* a 



dS v 
dy 

dy 



8m* 2 a 2 S. 



Am* 2 a 2 



qy — 



2F 

- ^-(35) 



The bulk solution of this equation is S* v = and 
Sl v = Sb — 2reEN F a, that coincides with the polar- 
ization obtained from (|26 p - (|2~7| at q — > 0. 

In order to calculate the spin polarization near the in- 
terface we employ the hard wall boundary conditions'^ 
for (|35p . Such boundary conditions can be easily ob- 
tained from Eq. (|3ip by performing its summation over 
k and integrating from y — to some point ya, placed at 
the distance much larger than I, but still small compare 
to l so . A simple analysis of Eq. (f3Tj) shows that up to 
o(a 2 ) the sum over k of the vector product in the l.h.s. 
of (|3"Tj) can be neglected, while the r.h.s and the term 
containing the electric field turn to zero identically. As a 
result, we get 



k y \y=Q ■ 



(36) 



According to ([3T| . the above sum is zero at y = 0. Hence, 
it is also zero at y = yo. The latter sum coincides with the 
spin current within its conventional definition 2 -, where a 
contribution associated with the charge density due to 
the second term of the velocity operator ([3]) is ignored in 
an electroneutral system. Using the gradient expansion 
of (|3~lj) this current can easily be expressed 25 through 
Siv\y=0i its y derivative and the last term in the l.h.s. 
of pip . In this way one arrives to the boundary con- 
ditions from Ref. 6;.819|. We generalize these conditions 
by adding possible effects of the surface spin relaxation 
(see also Ref. fToh - These additional terms are character- 
ized by the two phcnomcnological parameters p y and p z . 
Finally we obtain 

D^^-\ v=0 + 2Dm*a(S» v (Q)-S b ) = -pM0) 



-D 



dy 

dsiM 

dy 



lw=o - 2Dm*aS* v (0) = -p y Sl(0) . (37) 



One can easily see from (l3"5j) , (j3"7) and (|3"T|) that at p x / y — 
the homogeneous bulk solution S* v = 0, S^ v (0) = Sb 
turns out to be the solution of the diffusion equation ev- 
erywhere at y > 0. Therefore, in this particular case the 
z-components of spin clouds from many impurities com- 
pletely cancel each other and there is no spin accumula- 
tion near the interface, in agreement with Refs. 
At the same time, when pi ^ the out of plane compo- 
nent of the spin density is not zero. In the case of weak 
surface relaxation pi <§; D/l so one obtains from J37p - (l3"5l) 



Sf v (0) = OMpyTeE 



1 



2nhD 



(38) 
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where we inserted h to restore conventional units. It is 
interesting to note that in such a regime of small enough 
pi the surface polarization does not depend on the spin- 
orbit constant. 



IV. SPIN-HALL RESISTANCE AND ENERGY 
DISSIPATION 

As it was defined in the Introduction, the interface 
spin-Hall resistance is given by 



R 



sH 



N F J 



(39) 



where j is the DC current density, j = aE 1 with the 
Drude conductivity a — ne 2 T/m*. The so defined spin- 
Hall resistance is closely related to the additional energy 
dissipation which takes place due to spin accumulation 
and relaxation near interfaces of a sample. Indeed, as was 
shown in Ref. [6j, the spin accumulation is associated with 
a correction to the electric conductivity of DC current 
flowing in the x direction. For Rashba SOI the correction 
to the current density has the form 



Aj(y) = - 



a 2 kl dS* 



4m* T 2 dy 



(40) 



This expression is finite within the distance ~ l so from 
the interface. After integration over y one obtains a cor- 
rection to the electric current 



induces in 2DEG a nonequilibrium spin density around 
a spin independent isotropic elastic scatterer. The z- 
component of this density has a shape of a dipole di- 
rected perpendicular to the external electric field, while 
the parallel to 2DEG polarization is isotropic. Due to the 
D'yakonov-Perel' spin relaxation, the spin density decays 
exponentially at a distance larger than the spin-orbit pre- 
cession length. Noteworthy, that such a cloud exists even 
in the case of the Rashba spin-orbit interaction when the 
macroscopic spin current is absent. We also calculated 
the macroscopic spin density near an interface from tak- 
ing the sum of clouds due to many scatterers and inde- 
pendently averaging over their positions. Surprisingly, 
in the case of the hard wall boundary, the so calculated 
spin polarization exactly coincides with that found from 
the drift diffusion or Boltzmann equations&i^. In this 
case the out of plane component of the spin polarization 
is zero, while the parallel polarization is a constant de- 
termined by the electric spin orientation^. Besides the 
hard wall boundary we also considered a more general 
boundary condition containing the interface spin relax- 
ation, or the spin leaking term. For such a general case 
S* 0. This polarization can be associated with the 
local magnetic potential, because the system attains its 
local equilibrium within the S z spatial variation scale, 
which is much larger than I. The magnetic potential, in 
its turn, is related to the DC electric current density via 
the interface spin-Hall resistance. The latter was shown 
to determine the additional energy dissipation due to re- 
laxation of the spin polarization near the interface. 



4m* T 2 



: ^.(o) ■ 



(41) 



The corresponding interface energy dissipation (per the 
unit of the interface length) can be expressed from (139)) 
and (HB as 



AW = AIE 



-a rR sH j 



In its turn, R s jj can be determined from Eq. ([38 
can be easily seen that AW > if p y > 0. 



(42) 
. It 



Besides conventional semiconductor quantum wells, 
the results of this work can be applied to metal adsor- 
bate systems with strong Rashba type spin splitting in 
the surface states^. In this case the spin cloud can be 
measured by STM with a magnetic tip. 



V. RESULTS AND DISCUSSION 

Summarizing the above results, within the drift diffu- 
sion theory we found out that the intrinsic spin-Hall effect 
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